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I. INTRODUCTION 



Spatial optical correlations, both quantum and classical, represent a fundamental resource for developing technolo- 
gies as quantum imaging [HQ, super-resolution Q, etc., which could open unprecedented opportunities in the field of 
metrology, positioning and imaging. Various protocols have been proposed [illisj and experimentally realized [T^ - |25j . 
Among them a great interest has been attracted by the so called ghost imaging [2^ (GI), based on the correlation 
\ in spatial intensity fluctuations (speckles [13 )• In this technique a light beam crosses (or it is reflected) by an object 
Qli< to be imaged. However, the beam that crossed the object is detected by a detector without any spatial resolution 
(bucket detector). The image of the object is retrieved when the bucket detector signal is correlated with the signal of a 
. spatially resolving detector measuring a light beam whose noise is spatially correlated to the previous beam (reference 



beam). The first demonstration of this technique was achieved with nonclassical states of light [28|, known as twin 
beams, produced by parametric down conversion (PDC). Then it was shown, both theoretically and experimentally, 
that this result can be achieved also with beam-split thermal light [2^ - [33 |. although with a smaller visibility. 

In practice this technique can be useful in the presence of phase distortions (for example when the beam crosses 
a diffusive medium, like fog), where intensity correlations with a second beam allow one to retrieve the spatial 
information [sij . 

QQ , Due to the conceptual and practical interest attracted to GI, a lot of works have been devoted to clearly theoretically 
' describe [111, Issl - fSTj and to improve this protocol fssl - li^ . The main parameters of GI discussed in these works are 
typical for any kind of imaging: the signal-to-noise ratio (SNR), which shows how well the image of the object is 
■ distinguishable from the background, and the resolution. The latter is related to the number of speckles (spatial modes) 
[ contained in the image and is basically the number of elementary details of the object that can be reconstructed in the 
• ghost image. Note that various terms are used in the literature: SNR [s^jlH, is sometimes called contrast-to- noise 
T— I ' ratio (CNR) [i^; number of spatial modes [13,113 is sometimes called the number of details of the image [ioj . 
^ . Despite the huge amount of literature devoted to GI, several points are still not clear concerning this technique. 
• ^ For instance, in all existing experimental works, GI is performed by measuring either normalized intensity correlation 
, functions (CFs) [iO,!!^ or intensity CFs with the background subtracted (also called CFs of intensity fluctuations) (sol - 
^ ' [12 . At the same time, up to recently, most theoretical papers considered only the simplest CFs, without normalization 
- - -' or background subtraction. In a recent theoretical paper, Chan et al. [46| for the first time considered both the previous 
procedures and showed that they lead to almost the same SNR value, while GI through the measurement of usual CFs 
has much smaller SNR. However, they only considered GI with bright thermal light and presented only the results of 
numerical simulations. In [1^ Erkmen et al. developed a theory of the SNR of GI in the analog detection scheme 
using the not-normalized intensity correlation, providing results when ac-coupling of the photocurrents is performed 
by introducing signal frequency filtering. 

In this paper, for providing the reader with a general discussion of second-order [sij GI in view of practical 
applications, we generalize the theory of Ref [46| by including into consideration, beyond thermal light ghost imaging 
(ThGI), the case of twin beams (TwGI) and assuming an arbitrary brightness of the light source. Moreover, we 
consider another protocol for GI that exploits the variance of the difference signal from the object and reference 



detectors [37[. It turns out that this method of GI provides the same SNR as the background subtraction and 
normalization of the CF for middle and high brightness of the source. 

We analyze the influence of important aspects like the brightness of the source, the losses, the number of spatiotcm- 
poral modes collected by the detector, and the resolution on the performance of GI in terms of SNR. 

Finally, in this work we compare the theory with the experiment, for two particularly interesting cases: mesoscopic 
twin beams (photon numbers per spatiotcmporal mode about 0.2) and bright pseudo-thermal light, both detected by 
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CCD cameras. The former is the best GI option if high intensities should be avoided (for instance, for the sake of 
not damaging the imaged object), and the latter is the simplest way to achieve high SNR if using bright light is not 
a problem. 

II. THEORY 

In general, in order to obtain a ghost image, a spatially incoherent beam interacts with the object and then is 
collected by a bucket detector without any spatial resolution, while the correlated beam, which does not interact with 
the object, is registered by a spatially resolving detector, namely an array of pixels. This procedure is repeated /C 
times. The image of the object is described by a parameter S{xj), where Xj represents the position of the pixel j in 
the reference region. Usually, S has the form 

S{xj) = f{E[Ni],E[N2{x,)],E[N^N2{x)],E[nl],E[N^{xj)],...), 

i.e., a function / that involves the correlation function i?[N^Af|(a;j)], (p, q > 0), of the observable corresponding to the 
total number of photons collected at the bucket detector and at the j-th pixel of the reference arm, respectively. Here, 
E[X] = represents the average over the set of K. realizations. For the sake of simplicity we consider an 

imaged object defined by two levels of transmission, T = 1 and T = 0. 

The signal-to- noise ratio of a ghost imaging protocol can be defined as the ratio of the mean "contrast" to its 
standard deviation (mean fluctuation): 

SNRs = I {Srn-So.,t) \_ 

where Sin and Sout are the intensity values of the reconstructed ghost image, when xj is either inside (T = 1) or 
outside (T — 0) of the object profile, respectively, and 6S = S — (S) is the fluctuation. The mean value, denoted 
by (.), represents the theoretical expectation value, which can be estimated in practice by performing space averages 
over the regions "in" and "out" of the ghost image. 

The ghost image can be reconstructed by exploiting different GI protocols, namely different parameters S. Although 
all the protocols exploit the correlations between the object beam and the reference one, not all off them give the 
same results in terms of SNR. In particular, in the following we analyze four different protocols based on using: 

i) the Glauber intensity correlation function, when S{x) = G'-^^x) = E[NiN2{x)]. 

Here, the quantum expectation value is obviously {S{x)) = {NiN2{x)), and the standard deviation can be evaluated 
directly as {5'^S{x)) = {d'^[NiN2{x)]) /JC. The photon numbers are considered as quantum operators and the quantum 
mean values are evaluated in the appendix, where we describe in detail our theoretical model. 

ii) the normalized intensity CF, when Six) = g^'^\x) = G'^'^^ / {E[Ni]E[N2{x)]). 

We note that, unlike all the other parameters considered at point i, iii and iv, the standard deviation of g*-^-* 
measurement was evaluated in our work by means of the uncertainty propagation of the quantities G^^\ E[Ni] and 
E[N2{x)]. 

iii) the covariance, or CF of intensity fluctuations, S{x) = Cov{x) = E[(Ni - E[Ni]){N2{x) - E[N2{x)])] = 
E[niN2ix)] ~ E[Ni]E[N2ix)]. 

In this case we have {Gov{x)) = !^^{{NiN2{x)) — {Ni){N2{x))). The coefficient depending on the number of 
realizations, namely the number of acquired images, is the usual one allowing the unbiased estimation of the theoretical 
correlation. If we assume that i?[Ni] and i?[N2] can be well represented by their expectation values (Ni) and {N2) in 
the calculation of the fluctuations (it happens ii JC ^ 1), we have 

(S^Covix)) ~ {6'^[5NiSN2{x)])/IC = {{[SNiSN2ix)]^) - {5NiSN2{x)f )/IC. 

iv) the variance of the intensity difference, 

S{x) = E[{Ni-N2{x)~E[Ni~N2{x)]f] = (S[N?]-S[Ni]2)+(£;[iV2(a;)2]-i;[7V2(a;)]2)-2(S[NiiV2(x)]-£;[Ni]£;[iV2(^)]). 

Whatever the parameter 5, it involves in some form the CFs, which are second-order intensity moments. Therefore, 
the evaluation of its fluctuations requires expressions for the first- to fourth-order moments of the intensity. These 
are calculated in the Appendix, where the details of our model are extensively described. 

We consider several parameters that influence the SNR. 

• The brightness of the source, which is related to the average number of photons fi per single spatio-temporal 
mode. 
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• The total number M of spatio-temporal modes collected by each element (pixel) of the spatially resolving 
detector in the reference arm. It is approximately given by the product of the number of spatial modes, 
AAsp = inax[ylpij;/ylco/i, 1], and the number of temporal modes, AAt — max [T^ei /7co/i 1 1] ■ Here, ~A.pix^ Tdet are 
the pixel size and the detection (integration) time, respectively, while Acoh, Tcoh represent the characteristic 
coherence area (roughly, the speckle size) and the coherence time of the source. From the viewpoint of the light 
statistics, the total number of modes is M > 1; it is considered equal to 1 when Apix < Acoh and Tdet < Tcoh- 
While through the paper we always consider a usual situation in which the integration time is much larger than 
the coherence time, the case of Apix < Acoh and Apix > Acoh sltc both analyzed in the text. We stress that M 
should not be confused with the number of realizations, for instance the number of acquired frames /C: M is 
the number of modes collected by a pixel in a single acquired frame, so that X = {N2) = ri2MiJ, represents the 
total number of photons detected in the pixel, further called the "illumination level" . 

• We take into account the overall transmission-coUcction-detcction efficiency of the two channels, < iy < 1 with 
j = 1,2, i.e. the probability to detect an emitted photon. According to our model, developed in the appendix, 
the situation in which the pixel of the reference detector is smaller than the coherence area {Apix ^ Acoh) 
corresponds to a reduction of the collection efficiency. Therefore, 772 includes also a factor that takes into 
account for the geometrical collection probability of a photon in a certain spatial mode, 772 = 772.0 * Tl2,coih where 
Tlcou ~ min[Apix/Acoh, !]■ 

• Finally, the number of spatial resolution cells of the reconstructed image is represented by i? = 
Ain/ T^3''X-\ApixTAcoh]i with Ain the area of the ghost image where T = \. When the pixel size is equal or 
smaller than the size of a single spatial mode, R is determined by the number of spatial modes, otherwise it is 
given by the number of pixels in the area Ain- The ideal condition for maximizing in principle the resolution and 
the collection efficiency is Apix — Acoh, that is, the pixel size should approximately coincide with the speckle 
size 0- 

As an example of a specific relevant case, in Tab. Uwe report the expressions for the SNR calculated for the ideal 
case of unity transmission/quantum efficiency and collection efficiency in both optical paths, 771 = 772 = 1. The SNR 
values are normalized to \/K,. 
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TABLE I: Expressions for the signal-to-noise ratio (SNR) of the reconstructed ghost image with thermal light (ThGI) and 
with twin beams (TwGI) for different protocols described in the text, in the lossless case (77 = 1). 

First we observe that both for TwGI and ThGI the performances of the scheme based on arc drastically worse 
than for the other methods because the SNR drops as SNRq2 1/R{R^ 1) and as 1/a/M (M ^ 1) with respect 
to the resolution parameter R and with respect to the number of modes detected by each pixel M . At the same 
time, all the other protocols (exploiting _g^, Cov and Var) scale as 1/\/R{R 3> 1) and asymptotically ~ const for 
M ^ 1. We report the theoretical dependence of SNRs on R in Fig. [6] and Fig 17] in Sec. IIII CI This fact makes 
the not-normalized and not-subtracted correlation function , albeit largely considered in the literature (see, for 
instance, Refs. (bsI [isl l46j). inappropriate for ghost imaging when compared with the other protocols. Thus, we do 
not investigate it further through the paper. 

Fig. [T] presents the SNR for all the protocols as a function of the illumination level when M = 1 and ?? = 1. In 
the typical situation of large R (R > 10) it turns out that SNRg2, SNRcov and SNRvar behave the same for large 
values of detected photons X ^ 1, approaching approximatively the upper value of (2i?)~^/^, although they reach this 
bound at different values of the brightness. In particular, for all ThGI protocols the bound is reached when X ^ 1. 
For TwGI based on the variance and covariance, the condition depends on the resolution, namely X ^ l/(2i?), while 
our calculation shows that g^"^^ reaches the flat region as soon as I ^ l/(2i?^). Concerning the performance of a 
very low-brightness source. Fig. [T] shows that the advantage of using twin beams is in general very pronounced but, 
even with the same source, some protocol seems to be more convenient than others. In particular, for ThGI all the 
protocols scale as cx T with the exception of the variance method, for which SNR approaches zero faster, i.e. cx 
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FIG. 1: Signal-to-noise ratio of the ghost image as a function of the illumination level X for different protocols investigated 
in the paper. The number of detected spatio-temporal modes in a single run and the resolution parameter are fixed to be, 
respectively, M — 1 and 7? = 100. 

For TwGI all the methods lead to the same asymptotic behavior cx X^^^. Therefore, we can conclude that for the ghost 
imaging of a complex object the three protocols exploiting SNRg2, SNRcov and SNRyar have the same performance 
at medium and high intensities (for the same source, quantum or classical), while for very low brightness, preference 
should be given to the normalized correlation function g^. 

Now we focus on the difference between the signal-to-noise ratio of TwGI, SNRtw, that is obtained in the limit 
of large R, and the corresponding SNRth obtained with thermal light in the same limit, by using the covariance 
as the parameter for the GI reconstruction. In order to perform a fair and probably more useful comparison in a 
practical situation, we already introduced the "illumination" X = rjMfi^ i.e. the total number of photons detected 
in the pixel in the single image. In Fig. [2lja) we report the quantum and classical SNR, for different values of the 
transmission/detection efficiency, balanced for the two channels, while fixing M = 1 (a single temporal mode and a 
single spatial mode are detected). Not surprisingly, the curve for the thermal light is not influenced by the losses, while 
the TwGI curve, in general showing a better performance than the classical one, for high losses approaches the ThGI. 
The advantage of twin-beam light is evident for a low-brightness source, i.e. when the number of photons per space- 
temporal mode /i ^ 1. On the contrary, when /i 3> 1, thermal and twin-beam states of light produce the same results 
(Fig. [2])- The physical reason can be understood with the help of Eqs. (jB)), ^ and ((8)) for the two-mode statistics. 
Both classical and quantum light have the same single-mode thermal fluctuations (S'^rij) ~ rjfi + rf' ji^ . However, 
the photon-number correlations for thermal light, {5ni5n2)TH = {ni){n2) = Tj^^J^, are relevant when the number of 
photons is large, {rij) = fi > 1, and do not include the shot-noise component oc fi, which remains uncorrelated. On 
the other hand, the correlations in twin beams are {SniSn2)TB = V^f^i^ + m)j which shows that even the shot-noise 
component of the fluctuations is correlated. Moreover, classical two-mode correlations, generated by a beam splitter 
obeying the Bernoulli statistics, are not reduced by the losses, which are again described by a beam splitter model 
(see the Appendix and the discussion after Eq. ([5])). 

It is also important to present (Figl2jb)) the unbalanced case, ryi > ?72, which in our model includes the typical 
situation in many experiments where the pixel is smaller than the coherence area. In this case, 772, co/; ~ -^pix / ^coh < 1, 
and the advantage of the TwGI is reduced. We observe that thermal ghost imaging is weakly influenced by the 
unbalancing. This indicates that, in order to make the best of the TwGI, in the photon-counting regime one should 
properly set the pixel dimension. 

In FigdJ the dependence of SNR on the number of modes M collected by a single pixel is shown. We see that 
the TwGI is insensitive to M while the SNR for ThGI is reduced for large M but can be recovered at higher values 
of illumination X. A discussion about that point can be found in the Appendix after Eq. pop . This means that 
ThGI performs best in the single-mode regime, while TwGI, even at low brightness, can be easily brought to a 
high-illumination level by accumulating a large number of temporal modes in a single frame, without any decrease 
in the SNR. However, provided that the source has the same parameters, i.e. the same brightness /i, the same M 
and therefore the same X, the advantages of quantum light can be of several orders of magnitude, even for large 
illumination X. 
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FIG. 2: Signal-to-noise ratio (SNR) of a ghost image as a function of the illumination level X for twin-beam and thermal-light 
GI. The number of detected spatio-temporal modes in a single run and the resolution parameter are fixed to be, respectively, 
M = 1 and R — 100. (a); the three dashed curves refer to different values of the balanced detection probability rji — r]2 ~ rj. 
(b): the unbalanced case where r;i = 1 and 7)2 = 0.9, 0.5, 0.1 from the higher to the lower curve, respectively. 



In conclusion, our discussion demonstrates that single- mode {M = 1) ThGI is the best solution, considering the less 
demanding experimental resources, when one is not limited in the brightness of light. On the other hand, in situations 
where a low light level, < 1, is needed, for instance because of the photosensitivity of the object, quantum light 
provides a larger SNR (even for a relatively large illumination level). 



III. THE EXPERIMENT 



For the experimental comparison of GI obtained with thermal light and twin beams, we used two setups. In the 
first case, GI was performed with bright pseudo-thermal radiation (more than 10^ photons per radiation mode). In 
the second case, nonclassical radiation was used, namely, squeezed vacuum obtained via type-II parametric down- 
conversion. 

In both experiments, we used the 'numerical mask' method suggested in Ref. [i^. We selected two regions, one in 
the object beam, where we virtually inserted the mask, and the other in the reference beam. In the bucket channel, 
we simulated the presence of a completely opaque mask, which fully transmitted light within a certain area Ain and 
fully absorbed light outside it. This was done by simply setting to zero the signals of all pixels outside of this area. 
This method, which is only possible with an array photodetector in the bucket channel, presents the advantage of 
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FIG. 3: Signal-to-noise ratio (SNR) of the ghost image as a function of the illumination level I. Here we consider an ideal 
lossless situation rji = ri2 — 1 and the resolution is R = 100. 

arbitrarily varying the shape and size of the mask without any technical difficulties and provides an ideal limit for 
this part of the set-up in a realistic situation. 

A. Experimental setup with pseudo-thermal light 

The experiment with pseudo-thermal light ghost imaging was performed with the setup shown in Fig. |4] j43l |. 
Second-harmonic radiation of a Nd:YAG laser with the wavelength 532 nm, pulse duration 10 ns, and repetition rate 
47 Hz was incident on an Arecchi's rotating ground-glass disc generating pseudo-thermal light [3]. The radiation 
scattered by the disc was split by a non-polarizing beamsplitter and then both output beams were registered by two 
different parts of the matrix of a commercial digital photographic camera Sigma SD14. The far-field speckle patterns, 
symmetrical in the two parts, were focused on the matrix by a lens with the focal length 50 mm. As a result, the 
typical speckle size on the matrix was about 0.230 mm. With the pixel size of the matrix being 7.8 /im, a single 
speckle occupied, on the average, an area of 30 x 30 pixels. This was confirmed by the measurement of the spatial 
intensity correlation function. 

The rotation rate of the disc was chosen in such a way that the speckle pattern did not change noticeably during 
a single pulse but changed completely from pulse to pulse. The camera was taking frames with exposure times 1/50 
s, which was much larger than the pulse duration but still less than the distance between the pulses. As a result, 
most of the frames contained the speckle pattern from a single pulse ([i^. Fig. |4|) and few were 'empty'. The total 
number of captured frames was 5000. In each frame, we selected square regions in the two symmetrical parts (Fig. 
6 of Ref. d^), further referred to as the object and reference channels. As before, in the object channel, the 'masks' 
were introduced numerically, by simply ignoring the readings of the pixels outside them. 

The intensity of the laser light was reduced by using a Glan prism; by measuring the signals from separate pixels 
of the camera versus the input intensity we verified that the camera was operating in a linear regime and was not 
saturated. It is important that the output data of the camera was written in the RAW format. The mean values of 
signals from separate pixels were on the order of 300, which indicates that the number of incident photons was at 
least as much as that. Hence, per single speckle, i.e. for a single radiation mode, we had photon numbers as high as 
3 • lO**. As the pulse duration was much smaller than the typical time of the speckle pattern fluctuation, and the pixel 
size was much smaller than the speckle size, the number of modes was M = 1. The transmission/quantum efficiency 
was quite low (less than 0.01), due to the beam attenuation, but, as discussed in section II, this is not important for 
ThGI. 
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FIG. 4: Experimental setup for realizing ghost imaging with pseudo-thermal light. 




FIG. 5: Scheme of the setup for realizing ghost imaging with twin beams produced by Type II parametric down conversion. 

B. Experimental setup with type-II PDC 

In contrast to most ghost imaging experiments reported in the literature that are realized with faint PDC, here we 
set an experiment in the mcsoscopic regime using a mesoscopic PDC source and a CCD camera as a photon number 
resolving detector. This represents an important step in view of practical applications of this method. The setup is 
reported in FiglH 

For generating the twin beams we used a type II BBO non-linear crystal {I ~ 7 mm) pumped by the third harmonic 
(355 nm) of a Q-switched Nd:YAG laser. The pulses had a duration of 7J, = 5 ns with a repetition rate of 10 Hz 
and a maximum energy, at the selected wavelength, of about 200 mJ. The spurious non-Gaussian components were 
eliminated by spatial filtering (a lens with a focal length of 50 cm and a diamond pin-hole, 250 fj.m of diameter) of 
the beam that was re-coUimated with a diameter of Wp = 1.25 mm at the crystal. After the BBO, the pump was 
stopped by two UV mirrors (with ~ 98% declared transmission at 710 nm), and by a low-pass frequency filter (~ 95% 
transmission at 710 nm). The orthogonally polarized down-converted signal and idler beams were separated by two 
polarizers (~ 97% transmission) and finally the far field was registered by a 1340X400 CCD (Princeton Pixis:400BR 
with pixel size of 20 /.im, quantum efficiency of 80% and readout noise 4 electrons/pixel). A mechanical shutter 
provided the synchronization with the laser so that each acquired frame corresponded to the PDC emission generated 
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by a single shot of the laser. The far field was observed at the focal plane of the lens with 10 cm focus in an 
' / — / ' optical configuration. Our setup was developed exploiting the potentiality of the quantum scheme [l^, which, 
according to the theory, should outperform the classical light scheme for the same brightness (number of photon per 
spatio-temporal mode). For optimizing the detection of quantum correlations, basically two requirements should be 
fulfilled. First, transmission/detection and collection efficiency must be large. We realized it by reducing the optical 
losses after the crystal and using a CCD camera with high quantum efficiency, while the collection efficiency can be 
kept high if the pixel area Apix is larger than the coherence area Acoh (approximatively a speckle), i.e. a pixel should 
collect more than one spatial mode |27| . The second point is that the detection noise must be low compared with the 
shot noise. For both these reasons, it is very efficient to perform hardware binning of the physical pixels. It consists 
of grouping the physical pixels in squared blocks, each of them being processed by the CCD electronics as single 
" superpixel" . The number of photons collected by a superpixel is the sum of photon numbers of each pixel, whereas 
the readout noise is just slightly increased with respect to the one of a single pixel. We chose its dimension comparable 
to Acoh, i-e. the pixel (hereinafter we discard the prefix "super") size was set to 240 x 240(^m)^. The "illumination" 
X, defined as the number of photons detected by a pixel in a single shot image (frame), was about I = 1900. The 
expected number of temporal modes Mt = Tp/Tcoh detected in one frame was 5 • 10"^, considering the coherence time 
Tcoh of PDC around one picosecond. The number of spatial modes detected by the pixel was AAsp = Apix/ Acoh — 4. 
A measure of the coherence area, basically determined by the pump transverse diameter, is the size of the spatial 
cross-correlation of the signal and idler intensity patterns that has been evaluated to Acoh ~ 120 x 120(/J,m)^. The 
total number of modes M = x M.sp turned out to be compatible with the level of excess noise £ = I/M ^ 0.12 
due to the thermal statistics in the single beam. Thus, the number of photons per single space-temporal modes was 
/i =X/{'qM) ~ 0.20. Under these conditions, wc recorded 4000 frames. 



C. Experimental Results 

After having discarded all the frames corrupted by incoming cosmic rays, as well as 'empty' frames in the ThGI 
setup, we started the analysis by defining two correlated regions in the object and in the reference arm. In the TwGI 
setup, the regions included 13 x 15 = 195 pixels, while in the ThGI setup, they contained 750 x 750 pixels. These 
regions were the same for all the frames. 

In the twin beams set-up, one of the main problems we had to cope with was a strong instability of the Q-switched 
laser power from pulse to pulse (about 14%). Since the expectation value of the number of photons per mode /i is 
proportional to the square of the power P of the pump (fj. oc sinh^ (const -\/P)), fluctuations of P lead to the fluctuations 
of fi. As a consequence, the temporal statistics on many pulses is characterized by the mean fi and the variance V{p). 
Indeed, using the equations presented in the Appendix, we have obtained the expectation value of Cov inside and 
outside of the mask taking into account the pump instability: 

Cov',^ ^ v'Mipil + + V{^l)il + RM)], (2) 



Cozi^t = 7fM'V{^l)R. (3) 

When M and R are large, as in our experiment, the fluctuations of the pump dominate and this introduces a 
nonzero covariance for pixels in the reference arm that are, in principle, uncorrelated with the mask. The overall 
result is a strong noise of "background correlation" that hides the PDC spatial correlation, essential in the ghost 
image reconstruction. As a matter of fact, this effect would lead to an increased noise on the ghost image and a 
consequent decreasing of SNR. 

In order to compensate for the instability of the pump, we normalized each frame by its average value taken over 
a certain CCD spatial region. This procedure completely removed the problem. 

In Fig|S] we present the obtained experimental values of the SNR versus the resolution parameter R of the mask 
for the case of TwGI. The data set is in agreement with the theory where the experimental parameters M , fi and ry 
have been estimated independently. However, the best fit is obtained when rj = 0.42 that is smaller than the value 
0.62 that was measured very accurately in [l6|. This discrepancy can be explained by the relatively small pixel size 
(reduced collection efficiency of correlated photons) and the presence of background noise generated by detection and 
straylight. A detailed analysis of such problems can be found in [l6| . We stress that in our experiment we exploited 
successfully the theoretical results pointing out that the multimode regime is the most appropriate for TwGI. In fact, 
by collecting M ^ Aigp- M.t — 4 • 5000 wc reached a level of illumination X = 1900 such that the experimental sources 
of noise (like electronic read-noise of the CCD and straylight ecc.) were almost negligible. For we only report the 
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Pixels Number (R) 

FIG. 6: Signal-to-noise ratio for TwGI as a function of the resolution parameter R. The three series of experimental data refer 
to the different correlation parameters considered in our analysis g^, Gov, Var. They completely overlap in this regime, as well 
as the theoretical curves (dashed line). The theoretical curve for G^{R) is also reported (dotted line). The SNR is normalized 
by the square root of the number of frames of the sample IC — 4000. 

theoretical prediction, which is close to zero. In accordance to this prediction, our experimental data involving 4000 
frames was insufficient to produce any ghost image with protocol. 

In summary, we have demonstrated ghost imaging with twin beams in a mesoscopic regime {fi = 0.2, M = 
Msp ■ Mt — 4 • 5000, 1 = 1900) by a photon-number resolving detector for such parameters of the source (number of 
photons per mode and the illumination level) for which thermal ghost imaging provides much worse SNR, according 
to the values reported in Fig. |3l 

On the other hand, in figure [7] shows the results of ThGI experiment. The illumination level was about I sa 300, 
photon number per mode ^ « 10^, the number of modes M = 1, and the resolution parameter R was varied from 
i? w 4 to i? « 200 by changing the size of rectangular 'numerical masks'. One can see that the experimental points 
agree well with the theoretical dependence for SNR, which, according to Table U is in this case similar for all three 
methods (based on variance, covariance, and .g^^-*) and given by SNR = \/2R. 

IV. CONCLUSIONS 

For providing the reader, in view of practical applications, with a general exhaustive description of ghost imaging, 
we have discussed both thermal and twin-beam GI for four different protocols, based on the measurement of the 
following values: the Glauber correlation function G'^^'' = (NiA'2), the normalized Glauber's correlation function 
^(2) = G(^V(Ni)(A^2), the covariance Cov = {SN1SN2), and the variance of the difference signal Var = ((5^(Ni — iV2)). 

A first significant result, in accordance with the analysis in (46j for thermal light, is that the G*-^^ method performs 
much worse than the others (the SNR falling rapidly with the resolution parameter, i.e. the number of speckles 
transmitted through the mask in the object channel, and with the number of spatio-temporal modes detected by the 
reference detector): thus, albeit often considered in the literature for the sake of simplicity, the G*^^-* method is not 
worth using in real applications. 

On the other hand, the other methods have similar performances. For all of them, twin-beam GI performs largely 
better than thermal-light GI at low illumination levels (few photons per pixel in a single shot), while they become 
equivalent at high illumination levels (many photons per pixel in a single shot). Furthermore, TwGI is insensitive to 
the number of modes collected by a pixel and can be brought to a large illumination level by accumulating a large 
number of temporal modes in a single shot, even maintaining low brightness, which means low photon number per 
spatio-temporal mode. This is what is realized in our experimental demonstration, in which about 20000 modes are 
accumulated in the 5 ns duration of the pump pulse. ThGI performs best in the single-mode regime. In summary, 
ThGI is preferable whenever one does not have limits on the brightness (e.g. photo-sensible samples), due to its 
simpler realization. 

Finally, we have demonstrated the perfect agreement of the developed general theory of GI with experiment in 
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FIG. 7: Signal-to-noise ratio for ThGI as a function of the resolution parameter R. The three series of experimental data refer 
to the different correlation parameters considered in our analysis g^, Cov, Var. They completely overlap in this regime, as 
well as the theoretical curves (dashed line). The theoretical curve for G^{R) is also reported (dotted line) but in experiment, 
no readable image could be obtained from our dataset. The SNR is normalized by the square root of the number of frames of 
the sample. 



some particularly interesting (and sometimes unexplored) regimes for both TwGI and ThGI. 



V. APPENDIX 



In this Appendix we calculate the correlation functions needed for evaluating the SNR of different GI schemes. 

Both for thermal and PDC light we first consider the single-mode statistics and then we extend it to a more general 
multi-mode case. Only later we introduce the bucket detector that integrates the light over all the spatial modes in 
the 'mask' ('object') beam. 

PDC (48l |. in the approximation of a plane and monochromatic pump beam of frequency ujp, can be described by 
an evolution operator U that is the product of independent operators, each acting on a couple of symmetric modes 
{q,n) and (-q, -^1) [ij: 



U = exp(-^/" H{t')dt'] =(^c^p\^a\iq,n)al{-q,~n)~h.c], (4) 
\ ^Jo J L - J 

where ifi is the detuning of the signal(l) and idler(2) photons with respect to the degeneracy frequency ujp/2, ±q 
are the symmetric and correlated transverse wavevectors, a| and are the creation operators of the signal and idler 
photons, respectively, and ^(q, fl) is the coupling coefficient. Therefore, it can be seen as a collection of independent 
processes where couples of symmetric modes are correlated and the output field operators for signal and idler are 

61 (q, n) ^ [/(C)ai(q, il) + V{04{^q, -17), 

fe2(q, = C/(0a2(q, + V{Oa\{-q, -17), (5) 

where U{^) = cosh(^) and V{^) = sinh(^) determine the gain of the process by their maximum value and the 
coherence time by their inverse width in fl. Here the phases are not considered because they do not play any role in 
the photon-number statistics. 

In our model we assume rjj (j = 1,2) to be the detection probabilities for the photons in the mode j. Here, rjj is 
a coefficient that includes transmission, detection (quantum efficiency, QE, of the detector) and collection efficiency 
Vj.coii of the photon emitted in a certain spatial modes. In particular, in our model the collection efficiency can be less 
than one when the pixel of the spatially resolving detector is smaller than the spatial mode at the detection plane. In 
that case, the probability to collect a photon emitted in a certain mode is ijj^coii = ■Apix/Acoh- All the losses can be 
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taken into account by a model in which a beam sphtter with transmission r]j is preceding an ideal detector (QE= 1) 
(49| . For the calculation of the photon- number correlation function nj — bjbj, this leads to a simple substitution 
bj — > ^Jrfjbj in the normally ordered expressions for the operators. The single-mode photon-number statistics as well 
as the two- mode correlation can be obtained through a straightforward calculation by the relations ([5]), with the input 
state being the vacuum one: 



(rij) rij^i, (6) 

("j> = (1 + + 6/i^?7f ) , 

{n^) = liT,, (1 + 14/i77, + 36A*2r?| + 2A^l^7^^) , 



(nin2) = fi{l + 2fi)rjirj2, (7) 

(n^n2) = ^77i7?2 (l + 6/^^7]i + ^(2 + 47yi)) , 

(n^n^) = ^7^17^2 [1 + 2V37?i772 + 6h^{t]i + r]2 + 47^772) + 2^(1 + 2?7i + 27^2 + 27/1772)] , 

{nln2) = mm [1 + 24^^/2 _^ 18^2^7^(1 + + 2/i(l + 67?i)] . 



For our purposes we focus on the correlation in the far field obtained as the focal plane of a thin lens of focal length 
/ in an / — / configuration and around degeneracy (51 = 0). Here, any transverse mode q is associated with a single 
position X in the detection (focal) plane according to the geometric transformation {2cf/ujp)ci — >■ x with c being 
the speed of light. We observe that the diffraction effects arising from the transverse finite dimension of the optical 
system (typically the pump profile in the twin beam generation process) lead to a nonzero transverse coherence area, 
Acohi in the far field, representing the uncertainty in the arrival point of the correlated photons. The fluctuation of 
the difference photon number is {6'^{ni — 77,2))th' = (S'^ni) + ((5^77,2) — 2{SniSn2) = 2ri^{l — 7/) for 7] = 7/1 = 772. Thus, 
in the ideal lossless ca se ( 7 7 = 1), a perfect correlation appears in the photon number since {5'^{ni — n2))TW = 
(sub-shot- noise regime |14l. Il5|'). 

In the case of thermal radiation we assume a set of independent modes (i=l,..) addressed to the input of a 
50% beam splitter performing the transformations bu = {ai + ayac,i)/V2 and 62,4 = (a^ — a^ad) / ^2, where a„ac is 
the mode at the vacuum input port of the beam splitter. The photon-number correlation of the beams at the two 
output ports are obtained by exploiting the rule ((a|^)^a^) = p\{a\aiY , for normally ordered moments of thermal light. 
For the single-beam moments, taking into account the losses, we found again the expressions in Eqs. ([6]). This is in 
agreement with the fact that each beam from PDC, considered separately, manifests thermal statistics. On the other 
hand, the correlations of photon numbers between thermal beams 1 and 2 are 



(711712) = 2/^^77x7/2, (8) 

{nln2) = 2//^77i772 (1 + 3/^771) , 

{nlnl) = 2/x277i 772 [l + Un^ium + 3/i(7/i + 7/2)] , 

(71^712) = 2/^77 1 772 [1 + 9771/7 + Un'^ril] . 



In contrast to the quantum case, the correlation between the photon numbers in the coupled modes is not perfect; 
hence, for thermal light we have — 77,2))t/i = 277/i, which represents the shot noise limit. At the same time, the 

excess noise ^ /x^ in the second equation in (jS]) disappears regardless of the losses. 

According to our experimental realization (but without loss of generality), we start by considering a situation in 
which both reference and object beams are detected by a spatially resolving array of detectors (the pixels of a CCD 
camera in our case). Each pixel collects a number M > 1 of independent spatio-temporal modes in the far field and 
the total number of photons detected by the pixel is Nj = T,f£^nj^i, with j = 1, 2. It is easy to show that the moments 
of the N distribution can be expressed in terms of single-mode moments in Eqs. ([5]) as 
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{N,) = M{n,), (9) 

(N^) = M(n2) + Af(M-l)(n,)^ 

(Nf) = A/(n3) + 3M(^/-l)(?^2)(nJ)+M(M-l)(A^-2)(n,)^ 

(AT/) = M + M{M - 1) (3 (n^)" + 4 (n^) (n,)) + 6M{M - 1){M - 2) (n^) {n^f + 

+ M{M - 1){M - 2){M - 3) {ujf . 

Since the modes of signal and idler fields are pairwise correlated, the correlators of the count numbers in two sym- 
metrical pixels (collecting Af correlated modes) are 



{N1N2) = M{mn2)+M{M^l){m){n2), (10) 
{NIN2) = M {nln2) + M{M ~l){{nl){n2) + 2{nin2){ni)) 

+ M{M-l){M-2){nif{n2), 
{NlNi) ^ M{nlnl) 

+ M{M - 1) [{nl) {nl) + 2 {n\n2) (na) + 2 (yims)' + 2 {nml) (m)) 

+ M{M - \){M - 2) [(nl) {n2f + {nif (nl) + 4 (nms) (712)) 

+ A/(A/ - 1)(A// - 2)(A./ - 3) {nif (na)' , 

+ A/(A/ - 1) {{nl) (n2> + 3 (nj) (nms) + 3 (n^ij) (^i)) 
+ M{M - 1)(A// - 2) (3 (nin2) (ni)^ + 3 (nj) (m) (712)) 
+ M(Af - 1){M - 2)(A/ - 3) {jiif (n2> • 

The correlators (iViiV|) and (A^iA^I ) can be obtained by simply exchanging the indexes 1 and 2 in the second and 
fourth equations of pUj) . From the first lines of Eqs. ^ and (|10p we see that {SN1SN2) ~ M{SniSn2) (771 = 772 = 1)- 
By substituting the two-mode correlators following from Eqs. ^ and ^ we have {SNiSN2)tw = (^1) + {Ni)/M 
and {SNiSN2)th = {^D/M for quantum and thermal light, respectively. These expressions show that thermal 
correlations are sensitive to the numbers of modes in both channels, while the shot noise term is not reduced by the 
increase in M . This provides an intuitive explanation of the different behavior between TwGI and ThGI with respect 
to the number of modes M. 

Now we introduce the bucket detector whose output is represented by the operator Ni = Yj^^^^Nk.ii i-e. a sum of 
the photon numbers over R independent sets of spatial cells, each one collecting Af > 1 modes. R coincides with the 
number of pixels of the ghost image if the pixel is sufficiently large to collect more than one spatial mode, otherwise 
it coincides with the number of spatial modes (speckles) itself. 



(Nl) 


= R{Ni), 






= R{Ni)+R{R-l){Ni)\ 




{nl) 


= R (N^) + 3R{R - 1) (N^) (A^i) 


+ R{R - 1){R - 2) {Ni)\ 




= R{N^)+R{R-1){3{N^)^ ^ 


- 4 (Nl) (A^i)) + 6i?(i? - \){R - 2) (Nl) {N^f 




+ R{R-l){R-2){R-3){Ni)'^ . 





The correlators between the readings of the bucket detector and the ones of an arbitrary pixel in the reference channel 
registering light correlated with the one passing through the mask (dubbed with the subscript " in" ) are 
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{^lN2,^n) = {NiN2)+{R-l){Ni){N2), (12) 

(N?iV2.„) - {N^N2) + {R-l){{Ni){N2)+2{NiN2){Ni)) 

+ iR-l)iR-2){N,f{N2), 

- (mNi) + (i? - 1) (iVi) (N^) , 

+ (R-l) {(N^) (Ni) + 2 (iViiVl) (TVi)) 

+ (i?-l)(i?-2)(7Vi)^(7V|), 

+ (i? - 1) ((iVf ) (A^2> + 3 (A^f ) {NiN2) + 3 (iVi^A^s) (iVi)) 

+ {R-1){R- 2) (3 (7ViiV2> {Nif + 3 (7V2) (iVi) {N2)) 

+ {R - 1){R - 2){R - S) {N^f {N2) , 

- (N.Nl) + {R-1) (iVi) (7V2>' . 

The analogous equation for the "out" case is obtained by using the statistical independence of Ni and N2^out, i-e. by 
substituting (Nf A^l^^^) = {K){^lout)^ and then applying Eqs. ([11]) and ©. 

By using the definition ([1]) and Eqs. ([TT]) . ([T^ we arrive to our full analytical expressions for the SNR for the 
protocols considered in the paper. 
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